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1. Introduction 

Rosenfeld's fundamental measure theory (FMT) is a singularity in the world of 
approximate density functional theories. While all approximate functionals are built 
aiming at incorporating as much information on the uniform phase as there is available 
(Evans 1992), FMT is constructed on purely geometrical arguments (Rosenfeld 1989). 
Because of this, typical density functional recipes provide simple functionals with great 
flexibility to incorporate data of very different nature on the thermodynamics and 
the structure of the fluids, while fundamental-measure (FM) functionals have a very 
rigid structure which rejects almost any deviation from orthodoxy (Tarazona 2002, 
Cuesta et al 2002). The counterpart of this is that FM functionals exhibit a set of 
very special properties not shared by any other approximate functional. To begin 
with, the structure of the fluids is predicted rather than input (as in the other 
theories). Furthermore, FMT is naturally formulated for multicomponent systems, 
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while other theories have serious difficulties to pass from one-component fluids to even 
binary mixtures (Denton and Ashcroft f99f , Choudhury et al 2002). But perhaps 
the most striking and characteristic property which distinguishes FM functionals is 
dimensional reduction. This means that if a <i-dimcnsional system is constrained 
to lie in a <i'-dimcnsional subset (with d 1 < d) and we evaluate the d-dimensional 
FM functional at the density profile describing this d'-dimensional confinement, then 
we obtain as a result the FM functional for the d'-dimensional system (Roscnfcld 
et al 1996, Rosenfeld et al 1997, Cuesta and Martinez-Raton 1997a, Cuesta and 
Martinez-Raton 1997b). This extraordinary consistency of the theory (as a matter 
of fact at the origin of its rigidity (Cuesta et al 2002)) cannot be found in any other 
density functional theory, and is obviously a feature that exact functionals possess. 
Dimensional reduction can be extended down to zero-dimensional (0D) confinements 
(cavities holding no more than a particle), and as theory has developed, it has become 
more and more clear that this latter fact can be reformulated as a constructing 
principle of any FM functional (Tarazona and Rosenfeld 1997, Cuesta and Martinez- 
Raton 1997a, Tarazona 2000). 

FMT was first formulated as a continuum theory, but in a series of recent 
works the authors have extended the theory to lattice models and formulated a 
lattice fundamental measure theory (LFMT) (Lafucnte and Cuesta 2002, Lafucnte 
and Cuesta 2003, Lafuente and Cuesta 2004, Lafuente 2004). The extension has 
been most revealing about the structure of FM functionals. In fact, LFMT, whose 
roots are very close to the continuum theory (Lafuente and Cuesta 2002), has its 
most natural formulation as a cluster theory, becoming a kind of density functional 
version of Kikuchi's cluster variation method (Kikuchi 1951) in Morita's formulation 
(Morita 1994). In its latest formulation, LFMT can be constructed, for hard- 
particle models, out of the exact free-energy functional of a given set of 0D cavities 
(Lafuente and Cuesta 2004, Lafuente 2004). The final result is of an extraordinary 
simplicity, given the fact that it provides the exact functional of many onc-dimcnsional 
systems, and a good approximation for higher-dimensional ones (typically a Bcthe 
approximation for lattice gases with nearest-neighbor exclusion). On the other hand, 
LFMT exhibits dimensional-reduction consistency even beyond dimensional reduction 
(that is, we remain within LFMT under transformations more general than lower- 
dimensional confinements of the system; see (Lafuente and Cuesta 2003) for some 
examples, as well as section of this article). 

In this work we will present a formulation of LFMT based on a powerful 
combinatorial tool known as the theory of Mobius functions (Rota 1964, Aigner 1979, 
Stanley 1999), whose natural context are incidence algebras defined on partially 
ordered sets (or poscts). This formalism will allow us to rigorously prove a list 
of results about LFMT. Thus, after revisiting LFTM as reported in (Lafuente and 
Cuesta 2004) and analyzing its structure (section EJ, we will prove (section that 
given a lattice model as well as a basic set of clusters W ma x (to be precisely defined 
later) of the lattice, there exists a unique linear combination of the free-energy density 
functionals on subcluster of W max which yields the exact free energy when evaluated 
at 0D density profiles (i.e. density profiles with support one of those clusters). A 
special choice of W max gives rise to LFMT (section^, but the cluster expansion has a 
wider range of applications (section |5J) . We study the low-density limit of the cluster 
expansion and redefine the clusters of LFMT (0D cavities) in such a way that for 
any model it is guaranteed that the zero-density limits of the pair and triplet direct 
correlation functions are exact (sectional), and suggest how to extend the definition in 
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order to systematically incorporate higher order direct correlation functions. Finally, 
the behaviour of the cluster expansion under the action of a mapping is analyzed, a 
consequence of which is the proof that LFMT is closed under dimensional reduction 
and other more general mappings (section E}. 

The basic result used in the above proofs is a theorem about the Mobius function 
which is stated and proven in an appendix. The appendix also contains an important 
result of the theory of Mobius functions (the cross-cut theorem) that yields some 
simplifications in the calculation of the Mobius functions of a given poset. We conclude 
this article with a discussion about some of the consequences of this reformulation of 
LFMT (the simplicity of its application being perhaps one of the most remarkable 
because of its practical consequences) as well as some open questions related to this 
theory. 

2. Lattice fundamental measure recipe reviewed: a multicomponent 
example 

In this section we will review the procedure recently proposed by the authors (Lafucnte 
and Cuesta 2004) to construct a FM functional for any hard-core lattice model. As 
in the latest versions for continuum models (Tarazona and Roscnfcld 1997, Tarazona 
2000), the constructive principle is based on the exact dimensional crossover to 0D 
cavities. In brief, the aim of this procedure is to build the simplest functional (under 
certain assumptions) which applied to 0D cavities, produces the exact result. 

The first hypothesis of the recipe, based on the exact functional for one- 
dimensional hard rods (Lafuente and Cuesta 2002) and the common pattern shared 
by all lattice FM functionals studied by the authors (Lafuente and Cuesta 2002, 
2003), is that the excess (over ideal) free-energy functional of an arbitrary hard-core 
multicomponent system in a lattice C has the form (in units of kT, the Boltzmann 
constant times the temperature) 

seckei i=i teC « (a) 

where 1 is a set of indices which label the different weighted densities nS k > (s) ; au are 
integer coefficients which depend on the specific model; ^o{v) = 77 + (1 — 77) ln(l — rj) is 
the excess free energy of a cavity admitting no more than one particle in which < rj < 
1 is the average occupancy; p(s) = (pi (s), . . . , p p (s)) denotes the vector of one-particle 
density functions of the different p species, and C^ k \s) = (c[ k '(s), . . . ,Cp fe '(s)) arc 
vectors formed by subsets of lattice nodes, i.e. C^ k \s) C C (i is a species subindex). 
Hereafter, C^ k \s) will be referred to as (multicomponent) cluster or cavity. The 
kind of cavities involved in l|2.1|) are 1-particle cavities because they are such that if 
a particle of species i occupies a node of C^ k \s), it excludes all nodes of Cj k \s) to 
particles of species j, for any j = 1, . . . ,p (including i)L 

Associated to cavities are 0D density profiles. If C = (C\ , . . . , C p ) denotes a 
cavity, a 0D density profile associated to C is a density vector, denoted p c (s) = 
(pCi ( s )i ■ ■ ■ j Pc p ( s )) 1 with support C (i.e. pc t has support Ci). Note that for such p c (s), 

J In previous works this kind of object has been referred to as "0D cavity", but as we will redefine 
this concept later on, this new nomenclature is preferred. 
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the corresponding weighted density X)f=i J2tec Pd (*) IS ^ ne avera g e occupancy of the 
lattice subset defined by \J^ =1 C^ k \s). 

Given a specific model, the functional in i|2.1|) will be completely determined by 
the coefficients {a^} and the family of subsets {C^ k \s)}. As it was shown in a previous 
work (Lafucnte and Cucsta 2004), these unknowns can be uniquely determined by just 
imposing that the approximate functional (|2.1[1 recovers the exact limit for any 0D 
density profile of any 1-particle cavity. This condition can be expressed as 

•^fm [Pc] = ®o(v) f° r an y 1-particle cavity C, (2-2) 

with r\ = Yh=i J2tec PCift) average occupancy of the cavity C. From this, it is 
easy to notice that if the functional [p] satisfies condition (|2.2|l for the 0D density 
profile p c {s), then this condition immediately holds for any 0D density profile p c /(s) 
such that C' C C (inclusion here must be understood componentwise), simply because 
p C i is nothing but a particular choice of p c . Therefore, we can focus on the set of 
maximal cavities (Lafucnte and Cuesta 2002), which are cavities not contained in any 
other cavity. Thus condition (|2.2|) holds if and only if it holds for the set of maximal 
cavities, i.e. we can replace (|2.2(l by 

•^fmIPc] = ®o(v) f° r any maximal 1-particle cavity C. (2-3) 

In (Lafuente and Cuesta 2004) we showed (and this will be proved in full detail in 
this paper) that once we have determined the set of maximal cavities — which only 
depends on the geometry of the interaction — condition (|2.3|) completely determines 
the functional l|2.1|l for the given system. In other words, this condition uniquely fixes 
both sets {afc} and {C^ k \s)}. 

As an illustration of the procedure, we refer the reader to (Lafucnte and 
Cuesta 2004) for some specific examples of one-component models. Here, in order to 
review the recipe and to extend the collection of examples, we will apply the procedure 
to a one-dimensional binary hard-rod mixture. 

Let us consider a system in the one-dimensional lattice Z with two species of hard 
rods: the largest of length ox = 3 and the smallest erg = 2 (both in lattice spacing 
units). This model represents one example of the non-additive case which was exactly 
solved in (Lafucnte and Cucsta 2002). 

The first step of the procedure amounts to determining the set of maximal 
1-particle cavities of the system. As in previous examples, we will use a diagrammatic 
notation for cavities C = (Cl,Cs)- Remember that each Ci (i = L, S) is a subset of 
lattice nodes and thus we can associate Ci to the labeled graph whose vertices are the 
lattice nodes in Ci C C and whose edges are the bonds linking nearest neighbors in 
the embedding lattice C. In this representation we could have, for instance, 

C= ({s,s+l,s + 2},{s,s + l}) = (S-o-o,o-o) = CHD-o (seZ), (2.4) 

where we have merged the two components of C in the last diagram using different 
colors to denote each species (white for the large and black for the small) §. With this 
notation, the set of maximal 1-particle cavities for the hard-rod binary mixture is 

W max = { , S-o-o | s e Z}. (2.5) 

Since we want the functional in (|2.1|l to recover the exact limit for all 0D density 
profiles with support any cavity in W maX j there must be in (|2.1|) a contribution 

§ Notice that all nodes of the graph carry their corresponding label; the fact that we write a single 
label s on one of the nodes of the graph simply aims at avoiding clumsy notations. 
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for each C € W max . Notice that the weighted density n( fe '(s) can be identified 
with the associated cavity C^ k \s) = (c[\s), . . . ,Cp (s)) and the latter with the 
corresponding diagram. Therefore, we can write this contribution to i|2.1|) as 

y^[q(o-o-o)<l>o(c-c-o ) + a(o-c-c)$ (g-c-c )] ■ (2.6) 

sGZ 

Furthermore, if we want the functional to yield the exact OD limit, coefficients s((hhi ) 
and a(o-c-e) must be equal to 1 (this is always true for the coefficients associated 
to maximal cavities). This first analysis provides us with an initial guess for the 
functional of the system, namely 

^\p]=^2[M^) + M^>)]' ( 2J ) 

s6Z 

If this were the final functional, then it should satisfy condition (|2.3fl . If we take, for 
instance, p c , with C = e-c-o , as a test OD density profile and evaluate the functional 
Jpyp we obtain 

^fmi[p„4_o] =$o(^>) +$0(0-0^) + $o(«^) + *o(«-o) 

+$o(<^>) + $ (ch-o) +*o(«t-i) + *o(°t-i) +2$ ( Ot + 1 ), (2.8) 
where we have used that $o(0) = and that p c is zero outside C = o-c^ (which 
implies, for instance, that <E>o(«>-^-°) = ^o(«>-°))- Now, if we check <|2.8|) with 
condition l|2.3|l . we note that apart from the exact contribution $o(o-c-o) we also 
have a few spurious terms. Therefore the initial guess 1|2.7[) needs to be modified in 
order to eliminate these terms. Taking into account the general form (|2.1|) and the 
procedure explained in (Lafuente and Cuesta 2004), we propose as a second guess 

^fm 2 [p] =E[ $ o(^)+ $ o(^)+ fl ( <> ^) $ o( c ^ > )]' ( 2 - 9 ) 

sez 

where the coefficient a(o-e-o) has to be determined. The only way to remove the 
term $o(o-e-o) m (|2.8[) is to set a(o-f^o) = —1; thus 

•^FM 2 [P^U] =*o(«4-o) + $o(«^) + <^o{^) + M°t - l) +*o(o t + i). (2.10) 

We can now iterate the procedure to remove the term $0 ( o-e ) . The next guess is 

+ $o(S-e-o)-$o(o-c-o) + o(o-c)$o(e-e)]- (2.11) 

sez 

Choosing a(c-c) = — 1 we get -?fm3 [P c_J>-o ] = $o(o-e-°), and the exact limit is 
recovered for any OD density profile with support any maximal 1-particle cavity in 
the set { e-*^> I s e Z} C W max . By symmetry, it is easy to verify that the functional 
l|2.1ip needs no additional terms to fully satisfy (|2.3|) : therefore the functional 

^fm[p] =E[*o(^>^) + *o(^o) -*o(«^o) -*o («-•)] (2.12) 

sez 

is exact when evaluated at any OD density profile. Moreover, this functional coincides 
with the exact one, obtained by Lafuente and Cuesta (2002) through a different, more 
involved method. 

As it was discussed in (Lafuente and Cuesta 2004), there are two remarkable 
features in this procedure. First of all, if we start off the iteration with the terms 
associated to maximal 1-particlc cavities, then all clusters defining the weighted 
densities n^ k '(s) in (|2.1() are also 1-particle cavities of the system. Not only that, they 
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are intersections of the maximal cavities. Secondly, once we adopt the first ansatz 
•^FMiWi there is a unique functional which fulfills condition 12.311 (in other words, the 
specific scheme one follows in order to remove the spurious terms is irrelevant). This 
last statement will be rigorously proved in the next section. 

3. A cluster expansion for the free-energy density functional 

If we analyze the general expression (|2.1|) we notice that the FM excess free-energy 
functional is built from 1-particle cavity contributions. Furthermore, if we considered 
the system defined only in one of these cavities, say C^ k \s) = (c[ k \s), . . . ,Cp k \s)), 
the exact excess free energy would be ^o{ n< " k \ s )) ■ Therefore, if .F^p] denotes the 
exact excess functional of a system defined in cavity C, we can rewrite (|2.1|) as|| 

FfmIp}= E ^^cip]- (3-1) 

C cavity 

Using this cluster notation, the exact excess free-energy functional of the system 
under consideration can be denoted TJi\p], where C is the (multicomponcnt) cluster 
(£, . . . , C). Now taking into account that, in general, functional l|3.1|) is approximate, 
the exact one can be written 

mp]^mp}+ E «( c )^c x h (3.2) 

C cavity 

where ^ [p] is the error of the FM approximation (obviously an unknown functional) . 

As mentioned in the previous section, the coefficients a(C) are uniquely 
determined by condition (|2.3J) . which, as we have shown, can be implemented through 
an iterative procedure determined by the inclusion relations of the intersections of 
the cavities involved (considered as subsets of lattice nodes) . In this section, we will 
show that expression l|3.2f> is a particular case of Mobius inversion formula, a major 
result of the theory of partially ordered sets. The reader is strongly advised to consult 
the specialized literature on this subject (Rota 1964, Aigner 1979, Stanley 1999). 
Here, we will just introduce the necessary mathematical background to provide a 
comprehensible and rigorous foundation to LFMT. 

3.1. Mobius inversion formula in a nutshell 

A partially ordered set or poset is a set together with a partial order relation denoted 
< i.e. a binary relation satisfying rcflcxivity, antisymmetry and transitivity). An 
example of poset is the set of clusters involved in functional (|2.12|) : 

W = {c— t^o , o— o-e , o-e-o , c— c | S 6 Z}, 

with the order defined by C < C 1 if and only if C C C' . In general the order is 
partial because there may be pairs of elements which are not comparable. In our 
example neither e-c-o < o- e-o nor o-<x> < o-<x> arc true, while, for instance, 

s , s 

o— c— o <■ c— c— o . 

Let W be a poset. For any x, y € W such that x < y we define the interval 
[x,y] = {zeW\x<z< y}. A poset is said to be locally finite if all its intervals are 
finite. The set of all intervals of poset W will be denoted by Int(W). 

[| This rewriting of 12.11 is particularly interesting because the fact that any explicit reference to 
the function &o(v) nas disappeared makes it applicable to cavities more general than the 1-particle 
cavities considered so far. More on this later. 
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Given a locally finite poset W, the incidence algebra I(W, K) of W over field K is 
the set of all mappings /: Int(W) — > K (we will write f(x, y) for f([x, y])) with the usual 
vector space structure and the inner product (/ * g)(x, y) = Y^ze[x y] f( x i z )d( z i v) (/ 
and g being elements of 7(W,K)). Note that this product is well-defined because 
being W locally finite the sum contains a finite number of terms. It is easy to 
check that 7(W, K) is an associative algebra with (two-sided) identity S(x, y) = 1 
if [x, y] = [x, x] = {x} and otherwise. Another useful function of /(W, K) is the zeta 
function ((x,y) = 1 for all [x,y] G Int(W). 

An important result for incidence algebras is (the present statement is a simplified 
version of Proposition 3.6.2 on p. 114 of (Stanley 1999)): 

Proposition 1 Let / be an element of 7(W,K); then / has (two-sided) inverse (i.e. 
there exists f- 1 G I(W, K) such that / * f- 1 = f' 1 * f = 5) if and only if f(x, x) ^ 
for all x e W. 

It follows from this proposition that the zeta function £ of a locally finite poset 
W is invertible; its inverse is the Mobius function of W and is denoted by /^yv- This 
function can be obtained recursively from the definition of the inverse: 

Recursion 1 (C * £*vv = ^) 

Hw(x,x) = 1, for all ie G W, 

Mw(Z)#) = - X! ^w(z,y), for x < y with y e W. (3.3) 

x<js<y 

Recursion 2 * ^ = 5) 

/itw(x,a;) = 1, for all ir 6 W, 

Hw(x,y) = — /iyy(x,2:), for x <_ y with x,y G W. (3.4) 

In both Recursion 1 and 2. we have used a dot to indicate the summation variable. 

We are now ready to formulate the key theorem for the rigorous foundation of 
LFMT (from Stanley (1999), Proposition 3.7.1 on p. 116): 

Theorem 1 (Mobius inversion formula) Let W be a poset such that for every x G W 
the subset {y G W | y < x} is finite. Let /, g: W — > V(K), (V(K) being a vector space 
over field K). Then 

f{x)=Y,9{v) forallx-GW (3.5) 

y<x 

if and only if 

9(x) = J2f(y)^(y, x ) forallxGW. (3.6) 

y<x 

In the following, we will show that Mobius inversion formula can be applied to 
obtain a cluster expansion of the exact free-energy functional of a general lattice gas 
with arbitrary interaction. After that, it will be straightforward to prove that LFMT 
amounts to taking a particular truncation of this expansion. 
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3.2. Cluster expansion of the free-energy functional 

The cluster expansion we will propose here is based on the formulation of the cluster 
variation method by Morita (1994). There, Mobius inversion is used to approximate 
the entropy of a lattice model by a linear combination of the exact entropies of the 
same model restricted to a family of lattice clusters. This is exactly the same idea we 
find in approximation l|3.1|l , but here it is applied to the free-energy density functional. 

Let us consider a multicomponent lattice gas with underlying lattice C. Let 
us assume that we have a poset W whose elements are (multicomponent) clusters 
of lattice L (as in the previous example), such that cluster C is in W. Note that by 
definition C < C for all C G W. Now, let us consider the mapping f(C) = Tc\p\, where 
C G W and Tc [p] is the exact free-energy functional of the given model restricted to 
cluster C. Particularizing x = C in 13. 6[) and taking into account that fiy\>(C, C) = 1, 
we have for the exact free-energy functional of the system 

Tc[p] = * £ [p] + i-Vw(e,C)]Fc[p\, (3.7) 
c<c 

where ^rfp] = !?(£) is an unknown functional. Expression (|3.7|l provides a cluster 
expansion of the exact total free-energy density functional of an arbitrary system. In 
general, functional ^ c\p] cannot be computed exactly, but for suitable choices of the 
cluster set W some of its properties can be derived. 

Let us make a particular choice: let us assume that W consists of C as well as 
every non-empty intersection of the clusters of certain set W max . We will then show 
that ^^[p] vanishes for every density profile with support any cluster in W max . (Note 
that, by construction, the clusters in W max are the maximal elements of W — {£} 
with respect to the order relation.) 

Let T> G W max and p^ a density profile with support T>. If C G W then we have 

Fc[Pt>} = ^cnx>[Px>], 
with TzIpt^ = 0. By construction, if C n T> ^ then C n T> G W. Taking into 
account that !Fc[Pt>] = •^t>\Pt>\i fr° m l|3.7|) we obtain 

*c[Pt>] = ^w(C,£)^cn^[p c ] = J2 T£ \Pv\ J2 Mw(C,£). (3.8) 
cew e<t> cnx>=£ 

Now, we can use corollary [Jin | Appendix A"| to show that 

Y Mw(C, C)=0 for all T> G W max and £ < C. (3.9) 

CnT>=£ 

Therefore, the unknown functional ^ c fulfills the condition 

^c[p c ] = for all C G W - {£}, (3.10) 

since for every such C there exists at least one T> G W ma x such that C < T>, and so 
p c is a particular case of p-p. 

Now, if we approximate the exact free-energy functional in (|3.7|l by the truncation 

•Wp] = E [-»w(C,C)]Fe[p], (3.11) 
c<c 

the previous result 13.10f) guarantees that this approximation is exact for any density 
profile with support any cluster of W — {£}. 
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At this point, we have shown that given a specific lattice model there exists an 
approximation of the free-energy functional, given by (|3.11[) , which is exact when the 
system is restricted to certain set of cavities (namely, those in W — {£}). But, is it 
unique? Or in other words, how many free-energy functionals of the form 

^a PP [p] = a ( C )^c[p] (3.12) 
c<c 

are there which are exact in W — {£}? To end this section, we will prove that the only 
one with such property is H3.11|l . Before going into the technical details, notice that, 
as in the proof of l|3.10|l . it suffices that the above condition holds for the clusters in 
Wmax (the set of maximal elements of W — {£}). 

Let us suppose that we have an approximation like (|3.12[) such that it satisfies 

•^a PP [Px>] = Fc[Pt>] (3.13) 
for every cluster T> 6 W max . Since J~c\Pt>\ = ^cr\T>[PT>]i we have from (|3~T2|> and 
(EH51) 

•WPc] = Tt\Pt>\ = a(T>)F-n[ P v) + ]T T E [Pv\ J2 ( 3 ' 14 ) 

£<T> CnT>=£ 

This condition should be satisfied regardless the functional form of Tc [p] , therefore 
<|3.14[) can be cast in the recursion 

a(T>) = 1 for all V e W max , 

a(C) = for all S <T>. (3.15) 

CnT>=£ 

By defining a(C) = —1, the first of these two equations can be rewritten as 
Scnx>=x> a (^) = 0, so both equations can be gathered in the single one 

a ( c ) = for a11 V G w max and € < C. (3.16) 

CnT>=£ 

From H3.15fl it is clear that this recursion has a unique solution. On the other hand, 
H3.16fl shows that it is formally identical to l|3.9|l . Therefore, since it is linear, the 
solution must be a(C) = A/iw(C, C) for some constant A. Choosing C — C shows that 
A = —1; thus 

a{C) = -fi w (C, C) for all CeW. (3.17) 

In other words, the only functional of type yj.l'Jfy which is exact when restricted to 
the set of clusters W - {£} is EHip . 

Let us summarize the main results we have obtained in the following 

Theorem 2 Given an arbitrary lattice model, a certain set of clusters W max and the 
poset W formed by all non-empty intersections of elements of W max as well as the 
cluster £, then there exists a unique functional of type (|3.12|) which is exact when 
evaluated at density profiles p c with support any cluster C € W— {£}. This functional 
is given by 

^a PP [p]= Yl [-»w(C,C)]Fc[p], (3-18) 
cew-{£} 

where the integer coefficients ^>v(C, C) arc defined by either rccursion^^c. JjOJ) or 
recursion [3 (ec. I|3.4() ). 



Cluster density functional theory for lattice models 



10 



For some special cases, formula <|3.18ll is exact. One of them occurs if Tc\p\ is a 
local functional of p(s), i.e. 

»=i seCi 

With this particular choice for Tc [p] , the "error" functional l|3.7fl becomes 

cew i=i sec 

where hi(s) = J^cew PwiC, £)xc; an d XCi( s ) = 1 if s 6 Cj and otherwise 
is the indicator function of the set Cj. If we now define the cluster <7j(s) = 

(») 

(0, . . . , 0, o s , 0, . . . , 0), the function /ij(s) can be rewritten as 
h(s)= ^2 Hw(C,C), 

CntTi(s)=a-i(s) 

and because of corollary ^ of |A*ppcndix A"| /i,(s) = for all s <E £ and i = l,...,p. 
Hence *^[p] =0. 

4. Lattice fundamental measure theory revisited 

The Mobius formalism developed in section [3] endows us a very powerful alternative 
procedure to obtain the FM excess free-energy functional for any lattice system. 
The approximation l|3.18|l expresses the free-energy functional as a truncated cluster 
expansion of the form 1|3.1|) . Given that the ideal part of the functional is local in 
p(s) and, as we have just shown, for local functionals this cluster expansion is exact 
(regardless the choice of W max ), the approximation is actually made on the excess free- 
energy functional. If we now take for W ma x the set of maximal 1-particle cavities of a 
given lattice model, then theorem [5] tells us that the FM excess free-energy functional 
will be given by 

^fm[p]= E [-l*w(C,C)]* (nc\p]), (4.1) 
cew-{c} 

where nc[p] = Y^teC- Pi(t)- This expression is identical to (|2.ip . because in 

W — {£} there will be clusters of different shapes (labeled by k £ X in (|2.1[) ) and all 
their translates (labeled by s € C in (|2.1[l l. But now, thanks to theorem [21 we know 
that this is the only functional one can obtain with this particular choice for W ma x 
which is exact when evaluated at any 0D density profile with support any cluster in 
W-{C}. 

A remark is in order here. If C is a cluster of W — {£}, and C is a translation of 
C, then /zyy(C,£) = nw(C ,jC) because the Mobius function [iy\/(x,y) only depends 
on the interval [x,y] (c.f. its definition in recursion 1) and the translation operation 
is an obvious order-preserving isomorphism between [C,C] and [C',C]. This justifies 
why the coefficients at in jU} are independent of s G C. 

With this new formulation of LFMT it is a simple task to obtain the unknowns 
in expression (|2.1|l . As an illustration, let us re-derive functional l|2.12|l . The first step 
is to fix the appropriate cluster poset W, which contains all non-empty intersections 
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Figure 1. The Hassc diagram of a finite poset W is defined as the graph whose 
vertices are the elements of the poset and there is an edge between elements x 
and y if x < y and there is no z £ W such that 2 < .z < y. If x < y then y is 
drawn at higher level than x. This figure shows the Hassc diagram of the interval 
[o s , £]. Nodes connected through a descending path are ordered by transitivity. 
Because of this reason, a Hasse diagram is a very practical way of visualizing the 
order in a finite set. 



of maximal 1-particle cavities in the set W ma x described in lj2.5J) . as well as the cluster 
C. Thus, 

W = { c— 0— o , o— c— c , o — c — o , c— e , c— o , o— c , o— o ,c s ,o s |sG/3}U{/^}. (4-2) 

The second and last step is to compute the Mobius function ^w(C, C) for every 
C E W — {£} (note that Hw(C, C) = 1). This can be easily done by resorting 
to recursion ^ The natural iteration prescribed by this recursion is to start with 
maximal clusters (those in W max ) and then to follow a decreasing path. Notice that 
in order to obtain fj,yy(C,C) we only need the values of the Mobius function for the 
clusters in the interval (C,C\. To carry on this task, it is useful to draw the Hasse 
diagram of the corresponding interval (see figure since it explicitly shows the order 
structure of the latter. 

In our example, we should start with v-c-o and o-c-c (remember that the Mobius 
function does not depend on s, so what follows holds for any s £ £). As both are 
maximal clusters, then ( c-»^> ,£]=( o-c-c , C] = {£} (by definition, this is always 
true for maximal clusters). By applying recursion ^ we obtain 

Hw( c-e— o , jC) = /Uyv( o— c-e , C) = —fJ,y\;(jC, jC) = —1. 

In decreasing order, the next set of clusters involves o-e-o and c-c , and we have 

( C— C— O , jC] = { C — C — O , O— C— C , £}, f C-C , )C] = { C-C-K3 , C— C— C , jC.J. 

Thus nw(c^-(yo ,£.) = /Jw((y-c , C) = 1. Then we find c-o and o-c. By symmetry, 
/iyv(°-^,£) = l^w(°-t> I and since 

( c— o , £] = { c— c , o— c— o , c— c— o , c— c— o , o— c— c , £}, 
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we will have /zyy( o-o ,£) =0. Next we have o-o and c s . The corresponding intervals 
are, respectively, 

I o— o , jCJ = | c— o , o— c , c— c , o— c— o , o— o— o , c— c— o , c— ©— o , o— c— c , o— c— e , L-\ 
and 

/ —I f S S S S S S S 3 S 

I o s , L,\ — | c— o , o— c , c— e , c— e , o— c— o , c— c— o , c— c— o , o— c— c , o— c— c , *-}. 
Therefore /iyy ( o-o , C) = and pw(t> s ,C) = 0. Finally, for o s we have 

/ r S S S S S S S S 

y o s , L,\ = \ Q s , o-o , o-o , c— o , m , o-c , o-c , c— c , c— c , 

S S SS S S S S 

o— c— o , o-c— o , o— c-o , c— c-o , c-c— o , c-c— o , o-c— c , o— c— c , o— c— c , L,\ 

which leads to /ihi(o s ,£) = 0. Substituting these values of the Mobius function in 
the general expression l|4.1|l . we recover functional (|2.12|) . 

Although in our example we have shown that recursion^is enough to compute the 
Mobius function, there are many alternative (more efficient) techniques which exploit 
the order structure of the cluster poset associated to the given lattice model and make 
use of specialized results of the theory of posets (see section 3.8 of (Stanley 1999) and 
section IV. 3 of (Aigner 1979)). In |Appendix A| we provide a few examples of these 
tools. 



5. Extending lattice fundamental measure theory 

In the previous section we have achieved a reformulation of LFMT based on thcorcm[5| 
Actually, we have shown that LFMT is a particular case of the approximation proposed 
in that theorem, where we choose the cluster set W max as the set of maximal 1-particle 
cavities. We will see in the next section that, for hard-core models, this choice is an 
excellent balance between accuracy and simplicity of the approximate functional. But 
we want to stress that the cluster expansion of the free-energy functional (|3.7|) is more 
general: it applies not only to hard-core models, but is valid for any lattice model. 

In order to illustrate this, let us consider the Ising lattice gas in an arbitrary 
lattice C. The interaction potential for this system is such that each lattice node can 
be occupied at most by one particle and two particles interact with an energy J if 
they are placed at nodes which arc nearest neighbors. If wc applied theorem |21 to this 
system with W max the set of maximal 1-particle cavities (which for this model is just 
the set of all single lattice nodes), the approximate free-energy functional we would 
obtain would be that of the site-excluding ideal lattice gas, i.e. 

^appM = ]T [p(s) \np(s) + (1 - p{s)) ln(l - p{s))} . 
sec 

This choice for W max is clearly inappropriate for this system, since it ignores 
the interaction between nearest neighbors. In order to account for it wc should 
go beyond the standard LFMT and take 2-particlc cavities. Thus, W m ax = 
{all pairs of nearest neighbors} and W = {£} U {all pairs of nearest neighbors} U 
{all single nodes}. The Mobius function takes the values /xyv({s,t},£) = —1 and 
pw({s},£) = q(s) — 1 for every pair of nearest neighbors {s,t} and every single 
node {s} of C, q(s) being the coordination number at node s. From theorem [21 the 
approximate free-energy functional will be 

all n.n. {s.t} s££ 
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where 

F{s,t} [p] = P(s) In [p(s) - (s, t)] + p(t) In [p(t) - p^ (s, tj\ 

+ [1 - p(s) - p{t)\ ln[l - p{s) - p{t) + ( S , t)] , (5.2) 

p( 2 '(s, t) being the joint probability of finding two particles at nodes s and t. For this 
model it is not difficult to show that it can be eliminated in terms of p(s) and p(t) as 

f2V , 1 + C[p(s) + p(t)] yj{l + C[ P (s) + p(t)]} 2 4C(1 + C)p(s)p(t) 
P {) {s,t) = * — , 

with £ = exp(— J) — 1 (J is the interaction energy between nearest neighbors in kT 
units). Note that T{ s }[p] can be obtained from l|5.2|l by setting p(t) — and is 
just $o(p(s)). This approximation is exact when the lattice £ is a Bethe lattice^ 
(even with node-dependent coordination numbers) and it is equivalent to the Bethe 
approximation for any other lattice (Bowman and Levin 1982). 



6. Low-density limit of the approximate functional 

A fundamental ingredient of many approximate density functional theories is the 
available exact information about the low-density limit (Evans 1992; see also Cuesta 
et al 2002 for an analysis of this limit in the construction of the weighted density 
approximation and the FMT). In contrast, in our formulation of LFMT we have input 
the exact functionals for the system restricted to a certain set of clusters and used 
a combinatorial tool (Mobius inversion) to construct an "optimal" functional out of 
them. The latter is so restrictive that the only freedom we have in the resulting 
approximation l|3.18|) is limited to the choice of W m ax- In this section we will prove 
that with an appropriate choice of this set we can assure the correct behaviour of the 
approximate functional in the low-density limit up to, at least, third order. 

The second functional derivative of the excess part of the exact free-energy 
functional yields the exact pair direct correlation functional 

C (% S1 S2) _ ^tri 

%^ s *) S p i{8l )S Pj ( 82 ) 



fij(si, s 2 ) 



(p-0), (6.1) 



1 + f*k( s l, S 3 )pk(s 3 )fkj(S3, S 2 ) 

k S3(E£ 

where /,j(si,S2) = e~ < ^ >i ^ Sl ' S2 - ) — 1 is the Mayer function, (f>ij(sx, s 2 ) being the 
interaction potential (in kT units) between a particle of species i at node S\ and 
another of species j at node s 2 (the interaction potential is assumed pairwise). If we 
compute the pair direct correlation functional from the approximate functional l|3.18(l , 
we obtain 

c ipp, l3 ( s i^2) = [~Pw(C,C)]c ( i f(s 1 ,s 2 \C), 
cew-{c} 

where we have introduced 

^ (Sl ' S2|C) = "^(, 1 )^(, 2 )' 

K In particular, it is exact for the one-dimensional lattice Z. 



Cluster density functional theory for lattice models 



14 



the exact pair direct correlation functional for the system restricted to cluster C. Since 
c\j (si, s 2 |C) is exact, we have from ()6.1f) that 



c\f(s 1 ,s 2 \C) ~ /y(si,« 2 )xc 4 («i)x^(s2) (p -> 0), 

with xc( s ) the indicator function of C. Thus, the approximate pair direct correlation 
functional satisfies 

4p P ,ij( Sl ' S2 ) ~/y0i> s 2) 51 [-^w(c,£)]xc I (si)xc 3 (s 2 ) (p->o). 

C£W-{£} 

Taking into account that ^iw{C 1 C) = 1 and denoting try (si, s 2 ) = <Tj(si) U <Tj(s 2 ), 
we can rewrite the above expression in the more suitable form 

cS,,«(«i,«2)~/i 3 -(«i,a3)(l- E Mw(C,£)) (p-0). 

cry(ai,S2)<C 

Now, inasmuch as <ry(si,S2) < C is equivalent to C n er^ (si, s 2 ) = <Ty(si, s 2 ), it is a 
direct consequence of corollary ^ in | Appendix A"| that, when p — > 0, 

(2) f fij( s i, s 2) if <T.y (si, s 2 ) is contained in any cluster of yy max , 



C app,ij( S l> S 2) 



otherwise. 



The Mayer function has the same range of the interaction potential, when the latter 
is short ranged, /ij(si,s 2 ) vanishes if a particle of species i at Si does not interact 
with a particle of species j at s 2 . Therefore, if we want the approximate pair direct 
correlation functional to recover the exact low-density limit, it must happen that for 
any pair of nodes S\ and s 2 for which fij(s\, s 2 ) ^ the cluster tTij(si, s 2 ) is contained 
in at least one cluster of W max . 

Let us consider the particular case of hard-core interaction. For these systems, 
two particles interact if and only if they overlap. Accordingly, we can define a OD 
cavity as any cluster such that if we place two particles in any pair of its nodes (of 
any allowed species), they necessarily interact (in which case the corresponding Mayer 
function is non-zero). From this definition it should be clear that if we let W ma x be 
the set of maximal OD cavities, the density expansion of the approximate pair direct 
correlation functional will recover the exact zeroth order. 

Also, let us reconsider the Ising lattice gas. We showed in section [S] that the 
choice of W max as the set of maximal 1-particle cavities gives rise to a very poor 
approximation of the free-energy functional. Then, we took W max as the set of all 
pairs of nearest neighbors in order to account for the interaction. This set contains all 
maximal OD cavities according to the new definition, and so we are certain that the 
free-energy functional (|5.1|l gives a pair direct correlation functional with the exact 
zeroth order term in the density expansion. The new definition of OD cavity that we 
have just introduced is thus suitable for any kind of interaction (whether hard or soft, 
repulsive or attractive) in the sense that choosing W max as the set of all maximal OD 
cavities guarantees the correct low-density limit of c^ ' ^ (si, s 2 ). 

Now, let us look at higher order terms in the density expansion. If we take the 
third functional derivative of the excess part of the free-energy functional, we obtain 
the so-called triplet direct correlation functional 

Cy 3 fc( Sl ' S2 ' S3 ) = ~ X f U f U I \ ~ fij( s ^ s 2)fjk(s 2 , S 3 )fki(s 3 ,S 1 ) (p -> 0). 

J dpi{s 1 )dp j (S2)dpk(S3) 
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Computing c a J yfc(si, s 2 , S3) from i|3.18|) yields (reproducing the arguments given to 
obtain c^ Pjij ), when p -> 0, 



where iTyn^Si, S2, S3) = «Xi(si) U <Tj(s2) U <Tfc(s 3 ). Therefore, we are guaranteed to 
recover the exact low-density limit of the triplet direct correlation functional if every 
cluster cr ijk (s 1 ,s 2 , s 3 ) for which the product f i] {s 1 , s 2 )fjk(s2, s 3 )fkt(s 3 , «i) docs not 
vanish is contained in at least one cluster of W max . Note that because of (16. 1|) . if this 
holds then the density expansion of the pair direct correlation functional is exact up 
to first order in p. 

Again, the choice of W ma x as the set of maximal OD cavities (according to the new 
definition) is enough to assure the correct behaviour of ijfc(si, S2, S3) in the low- 
density limit. To see this, just note that the product fij(si, S2)fjk( s 2t s 3)fki( s 3> s i) 
is different from zero only if two particles of the corresponding species placed at any 
pair of nodes {si,S2,S3} interact; in other words, only if the nodes belong to the 
same maximal OD cavity. Therefore, the approximate free-energy functional 
with W max the set of maximal OD cavities recovers the exact density expansion of the 
pair and triplet direct correlation functional up to first and zeroth order, respectively. 

All this analysis can in principle be extended to higher-order direct correlation 
functions. In general, we will have that i± in (si, . . . , s n ) will recover the exact 
low-density limit provided the cluster cr^ i?i (si, . . . , s n ) = UT=i a n ( s ^ s contained in 
at least one cluster of Wmax for any combination of {(ii, Si), . . . , (i n , s n )} for which the 
exact low-density limit is different from zero. In practice, this is a very demanding task, 
because guaranteeing the correct low-density limit of c^ p {i in (si, . . . , s n ) amounts 
to consider all diagrams in the virial expansion of the nth order direct correlation 
function, and the number of them grows exponentially. Clearly, the current definition 
of OD cavity is insufficient to provide the correct behaviour beyond n = 3 (except in 
very particular cases, like some one-dimensional systems for which the approximation 
becomes exact), so if we want more terms in the low-density expansion we have to take 
bigger maximal OD cavities. Bigger cavities means more particles in one cavity and 
an increasingly higher difficulty to obtain the exact free-energy functional for a single 
cavity. So, as usual with expansions, although a systematic improvement is possible, 
going beyond the lowest terms may be too involved in practice. 

7. Dimensional reduction 

One of the most remarkable properties of FMT is dimensional crossover. This property 
means that by confining the particles of a d-dimensional system to lie in a (d — 1)- 
dimcnsional subset we obtain the (d — l)-dimensional FM functional out of the d- 
dimcnsional one. In the case of LFMT, a typical example would be to start with 
the lattice gas with nearest-neighbor exclusion in the simple cubic lattice in three 
dimensions and constrain the position of the particles to the nodes of one of the 
coordinate planes in order to obtain an effective system equivalent to the lattice 
gas with nearest-neighbor exclusion in the square lattice. In this example, the 
real dimension of the system is reduced from three to two, and in so doing, the 




Tijk (S 1 ,S 2 ,S 3 )<C 
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Figure 2. The system on the left corresponds to the lattice gas with first- and 
second-neighbor exclusion in the square lattice. If the position of the particles 
is confined to the black nodes of the lattice, the system behaves exactly like the 
lattice gas with nearest-neighbor exclusion in the square lattice (represented in 
the right figure). 



FM functional of the three-dimensional system is transformed into that of the two- 
dimensional one (Lafucnte and Cuesta 2003). But this is only an instance of a more 
general class of mappings between different models. As another example, the FM 
of the lattice gas with nearest-neighbor exclusion in the square lattice can also be 
obtained from the one of the lattice gas with first- and second-neighbor exclusion in 
the same lattice, as figure [5] illustrates (Lafuente and Cuesta 2003). 

The aim of this section is to define a general mapping between models (of which 
these two examples are particular cases) and to prove that the approximation (|3.18|l 
is "closed" with respect to this kind of mapping, i.e. that the FM functional of the 
original model becomes the one of the transformed model under the action of that 
mapping on the density profile. 

Let us assume that the original system has an underlying lattice C and its 
approximate free-energy functional is given from theorem [2 by 

•^app[p]= Yl [-Vw(C,£)]Mp] (7-1) 

C£W-{C} 

for a given choice of W max - If we restrict the position of the particles of species i 
to the embedded lattice £■ C C (lattice here is a general term which refers to any 
subset of C, finite or infinite), the approximate functional of the new effective system 
can be obtained by specializing functional (|T. l|t to a density profile with support 
C = (£[,..., C' p ), i.c 

Kp P [p] = ^ PP [Pc'} = E [-t*M c ,£)]Fcnc'[Pc'], 

C£W-{C\ 

where p c > (s) = p(s) if s € C and is otherwise. Since T<z [p] = 0, if W' denotes the 
set of all non-empty intersections of the clusters in W with the cluster C, the above 
expression becomes 

K pp [p} = E T c'[p] E [-l*w(C,C) + 6{C,Cj\, (7.2) 
CeW cnc'=C 

where we have introduced the identity of I(W,K), S(C,C), to compensate for the 
inclusion of /zyy(£, C) — 1 in the sum. We can now apply corollary ffl hi I Appendix K\ 
to evaluate the constrained sum on the r.h.s. of 1)7.2(1 . If C is contained in some 
cluster of W — {£}, then we recover the result of section i.c. T' p [p] = 3~c'\p\ 
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(in other words, the approximate functional of the effective system coincides with the 
exact one). On the contrary, if C is not contained in any cluster of W — {£}, then 

K V M= E [-^{C',C')]T C ,[ P ]. (7.3) 
c'ew-{C} 

Note that if yV max denotes the set of maximal clusters of W — {£■}, then the set of 
maximal clusters of VV' — {£■'} is the set W max of all non-empty intersections of the 
clusters of W ma x with C Consequently, the set W coincides with the set of all non- 
empty intersections of the clusters in W max , as well as the cluster C . In other words: 
the approximate functional (|7.3|l is the one that would have been obtained directly 
from theorem [5] if it had been applied to the effective system (the original system 
constrained to C) choosing W max as the set of maximal clusters. 

The previous result shows that LFMT is closed under dimensional reduction. 
Note that if we have the approximate free-energy functional prescribed by theorem 2 
for a given system and a certain W max , then for all systems that can be obtained 
from it through dimensional reduction, the resulting approximate functional for the 
lower-dimensional system coincides with the one deduced from theorem 2 with W max 
obtained from W max as above. In LFMT, the prescription for W max for a given 
system is to choose all maximal OD cavities. In order for LFMT to be closed under 
dimensional reduction, W max should coincide with the set of maximal OD cavities of 
the lower-dimensional system. A little reflection will convince the reader that this is 
indeed the case (and it is so as well for other kind of mappings, like the one described 
in figurc[2Jl. In general, if a system S is transformed by the action of a mapping T into 
another system T{S}, and if its set of maximal clusters W max (S') is transformed as 
described above into the set W max = T{W max (5')}, then the approximate functional 
<|3.18[) behaves consistently under T provided the prescription to choose W max (5') is 
such that W max (T{S}) = T{W max {S)}. 

8. Conclusions 

In this work we have provided a rigorous foundation of LFMT based on the formalism 
of Mobius inversion in posets of lattice clusters. The free-energy density functional 
is thus expressed in the form of a cluster expansion: given a set of basic clusters 
Wmax, there is a unique functional of the form (|3.18|) which is exact for every density 
profile with support any subcluster of W max . The cluster expansion (|3.18|l requires 
the exact expressions of the free-energy functional on the clusters of W max . The low- 
density limit of the functional (|3.18() dictates a definition of the clusters of W max 
— the OD cavities, or those clusters such that if there are two or more particles in 
them they necessarily interact — which guarantees the exact zero-density limit of the 
pair and triplet direct correlation functions. This redefinition subsumes the previous 
version of LFMT, valid for hard-core models (Lafucntc and Cuesta 2004), and extends 
it to include any lattice model with short-range interaction (the Ising lattice gas is 
an explicit example). The Mobius function formalism also allows us to analyze the 
behaviour of functional H3.18J1 under mappings between models, and a consequence of 
this analysis is the proof that LFMT behaves consistently under dimensional reduction 
(or confinements of the density into lower-dimensional sets of the lattice) . 

Going beyond LFMT — to account for higher order correlations, for instance — 
is, in principle, possible, but requires choices of W max containing larger clusters. But 
then the exact free-energy functional of these larger clusters is required, what may be 
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too involved. In some cases, though — e.g. in some one-dimensional models — , it is 
known that LFMT as such is exact, so no improvement is required in those cases. The 
reason behind this fact is not clear to us yet, and it is certainly a matter that deserves 
further thought. 

Another interesting point concerns the continuum limit. For some models (e.g. 
hard cubes in a d-dimcnsional cubic lattice) this limit is feasible and in fact the limit 
functional coincides with the FM functional of the continuum model (Lafuente 2004) . 
But there are important cases, like hard spheres, which are not easy to obtain as a limit 
of discrete lattice models. If this were possible, the result would be a functional which 
would recover the exact zero-density limit of the pair and triplet direct correlation 
functions, something that the best current FM functional for hard spheres does not 
accomplish (Tarazona and Rosenfeld 1997, Tarazona 2000, Cuesta et al 2002). From 
the insight provided by the present analysis of LFMT we dare to say that Tarazona's 
functional for hard spheres (Tarazona 2000) is the limit of the functionals l|3.18[) of a 
sequence of lattice models with an incomplete choice of W max (incomplete in the sense 
that some maximal 0D cavities are not in W ma x)- This is another problem certainly 
worth exploring. 
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Appendix A. Some technical results 

In order to make this work self-contained, in this appendix we will collect and prove 
some technical results of the theory of posets (and in particular of cluster posets) 
that we have used along this article. We will only prove those which are original 
contributions, while for the rest we simply address the reader to the specialized 
literature, where complete proofs can be found. 

In section 13 we have proved that the approximate functional (|3.18|) is exact when 
it is restricted to any cluster of the set W — {£}. Also, in sections and [7| we have 
analyzed the low-density limit and the dimensional crossover of the functional (|3.18(l , 
respectively. In all cases, the key point was to evaluate a constrained sum of the 
Mobius function fi\v(C,C) for those C £ W satisfying C (~l T> = £, T> and £ being 
clusters of the lattice C not necessarily in W. This relation can be expressed in the 
more general form o~ + C = £, where cr + is a mapping from the cluster posct W to 
another cluster poset which, in general, is different from W. For all instances in this 
work, a + is just the intersection with a given fixed cluster of the underlying lattice, 
and hence is an order-preserving map. 

The advantage of the above interpretation of the constraint over the clusters in 
W is that it allows to use an important result of the theory of posets which relates 
the Mobius functions of two posets, W and V, if both are connected through a pair 
(cr, cr + ) of order-preserving maps, where a is a Galois function, defined as 

Definition 1 Let V and W be posets. A mapping o: V — » W is called a Galois function 
if there exists a function er + : W — > V such that 

(i) a, a are order-preserving maps; 
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(ii) a + ax > x for all 16V and aa + z < z for all z G W. 

Given two poscts V and W, their Mobius functions are related via a Galois 
function in the precise way expressed in the following theorem (which is a version 
of Theorem 4.39 in p. 173 of Aigncr (1979)): 

Theorem 3 Let <r:V-*}Vbea Galois function. Then, for all x G V, y G W, 

/V+w^i a+ y) if x G <7 + W and y G <rV, 



z£W,cr+z=x 



otherwise. 



This theorem is exactly what we need to derive a result which can be applied 
directly to evaluate the constrained sums that appear along this article: 

Corollary 1 Let W be a cluster poset with underlying lattice C which contains all 
non-empty intersections of certain cluster poset W max as well as the cluster C Let 
C be a cluster of lattice £ not necessarily in W, and let us consider the poset W of 
all non-empty intersections of the clusters in W with C Then, for all C G W, 



X Mw(C,£) 



cew,cn£'=c' 




if £'G:X>for allX»G W - {£}, 
otherwise. 



Proof. To apply theorem 3, we first need to rewrite the constrained sum in the 
statement of the corollary in a more convenient form. Since ^ W', the only clusters 
in W contributing to the sum are those whose intersection with C is non-empty. Let 
W denote the poset of such clusters; then we have that /zyv(C, C) = /iyy(C, C) for all 
C G W, since [C, £] in W is identical to [C, C] in W. Therefore we have the identity 



E 

cew,cnz:'=c' cew,<r+c=c' 



where (7+ : W -> W is defined as tr+C = C n C. 

Now, in order to apply theorem |3| to (|A.1|I we need a pair (a,a + ) of order- 
preserving maps a: W' — ► W and cr + : W — > W, a being a Galois function. JLet us take 
for cr + the one introduced in l|A.l() and let us define a as <rC = inf{C <E W\C C C} 
for all C G W . It is a direct consequence of the definition that a and a + both are 
order-preserving maps. Moreover, for all C G W we have 

(70-+C = inf {T> G W | (7+C G X>} < C 

because a + C cC, and for all C G W', 

cr+o-C = <t + inf{X> eW|C'cD} = inf{X>' G W \ C G X>'} = C 

Therefore, a and cr + satisfy condition (ii) of definition^and thus a is a Galois function. 

At this point, we can apply theorem[3]to the r.h.s of ijA.lfl . Since for this case 
a + W = W', for all C G W we can write 

\ iiw>{C',£) if Ceo-W, 



cew.cnr^c' 



otherwise. 



The last step of the proof just amounts to showing that C G erW if and only if C 
is not contained in any cluster of W — {£}. Note that the latter is equivalent to 
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aC = £, therefore if it holds then C £ aW'. Now, let us assume that C £ aW' , then 
there exists a cluster C in W' such that aC = C. Since a is order- preserving and 
C < C , we have C = aC < aC ', where only the equality aC = C can hold, and 
the proof is complete. I 

To end this appendix we will bring about some results which simplify the 
calculation of the Mobius function of certain posets, and we will provide two 
applications related to the cluster posets involved in this work. First of all, we will 
give some definitions relative to a special type of posets: lattices (a mathematical 
notion not to be confused with physical lattices), since they appear in a natural way 
when we have to compute the Mobius function of a locally finite poset. 

A poset V is a lattice if for any x,y £ V, sup{x,y} and inf{x,y} are in V . One 
instance of a lattice is any interval [C, C] of a cluster poset W such as those involved 
in theorem|21 (which are closed under non-empty intersections) . A 0-element of a poset 
W, denoted 0, is an element satisfying < x for all x £ W. Dually, a 1-clcment of W, 
denoted 1, is an element satisfying x < 1 for all x £ W. Obviously, any finite lattice 
has a and a 1. A point of a finite lattice V is an element satisfying < x for which 
there is no element z £ V such that z < x. A copoint of a finite lattice V is an element 
satisfying x < 1 for which there is no element z £ V such that x < z. A subset M. of a 
finite lattice V is called a iower cross-cut if ^ and for all x £ with x ^ A4 
there is an element y £ M. with y < x. Dually, a subset is an upper cross-cut if 
1 ^ M. and for all 1 ^ x £ 7> with x ^ M. there is an element y £ M. with x < y. An 
example of lower (upper) cross-cut is the set of all points (copoints) of a given finite 
lattice. 

In our particular case, we have to compute the Mobius function /^yy(C, C), where 
W is a cluster poset which contains all non-empty intersections of the elements of 
certain cluster set W max as well as the cluster C, which is a 1-element. Let us 
consider the finite lattice Vc = [C,C] with = C and 1 = C. From recursion ^ 
it is straightforward that for all C £ W we have /xyy(C, C) = /ip c (0, 1). Therefore, in 
all cases we have to compute the Mobius function /i-p(0, 1) of certain finite lattice V. 
This fact makes the following theorem (from Aigner (1979), Theorem 4.42 on p. 175) 
very useful: 

Theorem 4 (Cross-cut theorem) Let V be a finite lattice, M a lower (upper) cross- 
cut and nk the number of sets A C M. with k elements such that sup^l = 1 (inf A = 0). 
Then 

M6,i) = ]T(-i) fc n fc . 

fe>0 

When we take in this theorem the lower (upper) cross-cut as the set of all points 
(copoints) of the lattice V, a direct consequence is the following corollary (from Stanley 
(1999), Corollary 3.9.5 on p. 126): 

Corollary 2 Let V be a finite lattice with point set Q and copoint set 1Z. Then 

MO,i) = o 

if 6 ^ inf TZ or 1 ^ sup Q. 

Note that this result renders the calculation of some values of the Mobius function of 
a given poset straightforward. 
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A practical application of corollary [21 is to easily compute some values of the 
Mobius function, say, in the working example of section 0] Note that it is a direct 
consequence of this corollary that (j,yy ( c-o , C) = /i>v(o-o , £) = A < w(o-o,.£) = 
Hw(v s ,C) = £*w(°s,£) == for any se£, since for all these clusters the supremum 
of the sets of points of the corresponding interval is different from C (a glance at 
figure n is enough to realize it). 

An important consequence of corollary concerns the form of the cluster 
expansion of the free-energy functional. Notice that in our formulation of LFMT 
we have worked with W defined as the cluster C and the non-empty intersections of 
the clusters in W m ax- We have shown that this choice is enough to ensure our main 
purpose, building an approximate free-energy functional which is exact in the clusters 
of W max . Having this idea in mind we could have started with the cluster poset V made 
of cluster C and all non-empty clusters contained in some cluster of W m ax (note that 
W C V). Now, if we compute /iy(C, C) for all C £ V, then the application of corollary|21 
implies that if C is in V but not in W then /iv(C, C) = 0, while nv(C, C) = jUyv(C, £) 
otherwise. In other words, if a cluster C is not the intersection of maximal clusters, 
then it does not contribute to the cluster expansion of the free-energy functional. So, 
we can constrain this cluster expansion to W without loss of generality, as we have 
indeed done. 
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